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ON THE LIMIT ZERO DISTRIBUTION OF TYPE I 
HERMITE FADE POLYNOMIALS 


NIKOLAY R. IKONOMOV, RALITZA K. KOVACHEVA, AND SERGEY P. SUETIN 

Abstract. In this paper are discussed the results of new numerical 
experiments on zero distribution of type I Hermite-Pade polynomials of 
order n — 200 for three different collections of three functions [1, /i, / 2 ]. 

These results are obtained by the authors numerically and do not match 
any of the theoretical results that were proven so far. We consider three 
simple cases of multivalued analytic functions fi and / 2 , with separated 
pairs of branch points belonging to the real line. In the first case both 
functions have two logarithmic branch points, in the second case they 
both have branch points of second order, and hnally, in the third case 
they both have branch points of third order. 

All three cases may be considered as representative of the asymptotic 
theory of Hermite-Pade polynomials. In the first two cases the numer¬ 
ical zero distribution of type I Hermite-Pade polynomials are similar 
to each other, despite the different kind of branching. But neither the 
logarithmic case, nor the square root case can be explained from the 
asymptotic point of view of the theory of type I Hermite-Pade polyno¬ 
mials. 

The numerical results of the current paper might be considered as a 
challenge for the community of all experts on Hermite-Pade polynomials 
theory. 
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1. Introduction 


In this article, we are concerned with the general problem of the limit 
zero distribution (LZD) of type I Hermite-Fade (HF) polynomials for the 
collection of m function [/o = 1, /i,..., fm-i], m ^ 3, where /i,..., fm-i 
are multivalued analytic functions, each one having a finite set of branch 
points on the Riemann sphere C; the reader is referred at first to the basic 
surveys [52], [69], and also [27], [49], [3]. Up to the end of the paper, suppose 
that the functions /o = 1, /i,..., fm-i are rationally independent (over the 
field C(z)). 

We restrict our attention to the case of the collection of three functions, 
in other words, to the case m = 3. Let us recall some basic notions from 
the HF polynomials theory (see [52], [69]). 

Suppose that two functions fi and /2 are given by convergent power 
series at the infinity point z = oo. For a fixed number n G N, the type I 
HF polynomials Qn,o,Qn,i,Qn,2 of degree ^ n, Qn,j = Qn,j{z-, fi, f2) G Fn, 
:= Cniz), j = 0,1, 2, not all Qn,j = 0, are defined in such a way that the 
following basic relation holds true: 



It is well known that such polynomials Qnj always exist and in the “generic 
case”^ they are unique up to a suitable normalization (see [52], [69], [49], 


[38]). 


To be more precise, we treat here, from numerical point of view, three 
different pairs of Markov-type functions: /i,/ 2 , 91,92 and hi,h 2 - 

Recall that a function / G =^(oo) is a Markov-type function if f{z) = 
Jl{z) + const, where ^ is the Cauchy transform of a positive Borel measure, 
with support supp // d M (a compact subset of the real line M), 



( 2 ) 


Up to the end of the paper, a support of a Markov-type function / = 
/2 -|- const means the support of the corresponding measure 9 ,. 

Let El := [ 01 , 02 ] = [—l,a], E 2 := [5i,f>2] = [— a, 1], where a G (0,1) is a 
real parameter. Hence both segments Ei and E 2 are overlapping, EiriE 2 = 
[—o, o] / 0. 


^In other terminology, “in common position”. 
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Case 1. Let 


/i(^) : = 

dx z — a 

/ = log , 

Z i El, 

(3) 

f2{z) ■■ = 

dx , z — 1 

/ = log , , 

J-a z — x z + a 

z ^ E 2 . 

(4) 


We take the main branch of the logarithmic function in (3) and (4), in the 
sense that 

z - a z-l 

log-, log-Rilogl = U, as z —>■ oo. 

z + 1 z + a 


Case 2. Let 

9i{z) 

92{z) 


z — a\ 
z + lj 


z - 1 
z + a 


1/2 

1/2 



a — X dx 
X + 1 X — z 


+ 1 ) 


1 — X dx 
X + a X — z 


z ^ El, 


z ^ E 2 . 


( 5 ) 

( 6 ) 


We take such a branch of the function, that gi{z),g 2 {z) —>■ 1 as z —>■ 00 , 
and under the square root function -yW we mean the “arithmetic square root 
function”, that is = x for x G IR+. 

Case 3. Let 

/ _ \ 1/3 

( 8 ) 


We take such a branch of the function, that hi{z), h2{z) —)> 1 as z —>■ 

00 , and, in what follows, under the cubic root function -yW we mean the 
“arithmetic cubic root function”, that is = x for x G M+. 


Clearly, all three functions fj,gj, hj, j = 1, 2,3 are Markov-type functions. 
But since a G (0,1), we have Eir\E 2 = [—a, a] / 0 and Ei / £' 2 . Hence, all 
three pairs /i,/ 2 , 91,92 and /ii,/i 2 are neither an Angelesco, nor a Nikishin 
system (see [52], [21], [49], [26], [27], [4]). Thus so far, not even one theorem 
of general type is known, which predicts the LZD of type I HP polynomials 
for the three collections of multivalued analytic functions [1, /i, / 2 ], [1, 91 , 92 ] 
and [l,/ii,/i 2 ], that were introduced above. Hence, the numerical results 
that are presented in this paper should be explained from theoretical point 
of view by the forthcoming papers. 

Since o is a parameter, we actually have three families of Markov-type 
functions: ^ = {fi = fi{z;a), f 2 = f 2 {z-,a), a G (0,1)}, ^ = {gi = 
giiz-,a), g 2 = g 2 {z]a), a G (0,1)}, and Jif = {hi = hi{z]a), h 2 = 
h 2 {z;a), a G (0,1)}. Thus, for a fixed n G N, we have three families of 
HP polynomials defined by the basic relation (1). To distinguish between 
Case 1, Case 2, and Case 3, we introduce new notations for the HP poly¬ 
nomials for the collections of three functions [l,gi,g 2 ] and [l,/ii,/i 2 ]. We 
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set Pn,o, Pn,i, Pn ,2 £ \ {0} for the HP polynomials corresponding to the 

collection [ 1 , 51 , 92 ], 


{Pnfl • 1 + Pn,l9l + Pn,2g2){z) — O 



z —>■ 00 , 


( 9 ) 


and Unfl, Un,i, Un ,2 £ IPn \ {0} for the HP polynomials corresponding to the 
collection [l,/ii,/i 2 ], 


{Unfi • 1 + Un,lhi + Un,2h2){z) — O 



z —oo. 


( 10 ) 


Notice that for each value of the parameter a £ (0,1) the supports of the 
two Markov-type functions fj and gj, j = 1,2, coincide with each other, 
but the types of the corresponding branch points are of different nature. 
It is worth noting that in each pair fi , /2 and 5i, 52 , the branch points 
are separated from each other (cf. [54], [2]). But since for each a G (0,1), 
El n E 2 / 0 , El E 2 , and E 2 Ei, the system /i ,/2 is neither an 
Angelesco, nor a Nikishin system. The same is true for the systems 51,52 
and hi, /i 2 . 

In all of the numerical experiments discussed in the present paper, we set 
n = 200 (see (1)). We also fix the set of values of parameter a as to be 
a = 0.2; 0.4; 0.625; 0.73; 0.8.2 

At first we set a = —0.1 (and keep n = 200). Clearly, that value of a 
is out of the range of (0,1). But the reason is that for a = —0.1 all three 
pairs fi, f 2 , 5 i, <?2 and hi, /12 form Angelesco systems, since Ei,E 2 C M and 
El n E 2 = 0 . On Fig. 9 and Fig. 11 zeros of HP polynomials Q 200 ,j and 
.^ 200 , 1 , j = 0,1,2, are plotted. The numerical distribution of the zeros of 
the HP polynomials are similar in both cases and do not depend on the 
different type of branching of the functions fi , /2 and 51,52 at the points 
z = ±l,±a. The reason is as follows. All these numerical results are in 
full agreement with the general theory of limit zero distribution (LZD) for 
HP polynomials in the Angelesco case (see [21], [26], [4]; in these papers the 
case of type H HP polynomials was considered, but it is well known that in 
the Angelesco’s case there is a direct connection between LZD of type I and 
type H HP polynomials). Since the sizes of the supports of Markov-type 
functions fi and /2 (and 51 , 52 as well) are equal, the phenomena of the 
so-called “pushing of the charge”^ is absent in the case a = —0.1. It is in 
full agreement also with the general theory (see [30], [1], [21]). From Fig. 9, 
Fig. 10, Fig. 11, Fig. 12 follows, that the zeros of HP polynomials (5200 ,1 ^-nd 
( 5200,2 (and p 200 ,i, -^ 200,2 as well) are located on the segments Ei and £' 2 , 
respectively. But the zeros of Q 200,0 (and £ 200,0 respectively) are located 
on the imaginary axis. Thus in some sense the zeros of (5200 ,0 “intend to 
separate” the zeros of (5200 ,1 and ( 5200,2 from each other. Again this result is 
in good agreement with the well known LZD of type I HP polynomials in the 
Angelesco case (see, e.g. [54]). To finish the analysis of the case a = —0.1, we 
emphasize that in this case all numerical results are in good agreement with 
the general theory of HP polynomials. Thus, it is reasonable to consider all 
numerical results for a = 0.2; 0.4; 0.625; 0.73; 0.8 as to be trustable. 


^For Case 3 we set a = 0.85 instead of 0.8. 

^The phenomena was discovered in 1981 by Gonchar and Rakhmanov in [21]. 
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We describe in short the numerical experiments. 

We set the parameter a to be 0.2; 0.4; 0.625; 0.73; 0.8, and for each value of 
a we find numerically the zeros of HP polynomials Q 200 ,j{^'i cl), P 200 ,j{z]a) 
and U 2 ooj{z] a), j = 0,1, 2, respectively. Thus we obtain a “series” of numer¬ 
ical experiments for the different values of the parameter a. By the reasons 
mentioned above, it is meaningful to consider all these numerical results as 
trustable, and hence we can analyze them from a theoretical point of view. 

Before laying out in details the discussion and the analysis of the results 
of the numerical experiments (see subsection §4.3), we recall basic general 
facts about Fade approximants'^ convergence theory (PA-theory). For this 
purpose, we provide in the next sections §2 and §3 a survey for the most 
important results of PA convergence theory, which is intended to make more 
clear the mathematical essence of the numerical results presented in the 
current paper. 


2. Classical Fade approximants (or J-fractions) convergence 

THEORY 

2.1. Stahl’s 5-curve and Stahl’s limit zero-pole distribution theo¬ 
rem in classical Fade approximants theory. The basic general results 
of PA-theory are of two types. The first type is concerned with the limit zero 
distribution of Fade polynomials for multivalued analytic functions with a 
finite set of branch points on the Riemann sphere C. Commonly, these re¬ 
sults are referred to as “weak asymptotics of PA”. The second type of results 
are devoted to the so-called “strong asymptotics of PA”. 

The problem of first type was originally stated in a formal way by J. Nut- 
tall (see [51] and the review [52]) at the beginning of 1970’s, and for the clas¬ 
sical (“one-point”) Fade approximants the problem was completely solved^ 
by H. Stahl [64]-[68] in the middle of 1980’s (see also [71]). 

The second type problem was also initiated by Nuttall (see [50], [52], [53], 
[55] and also [31, § 2, Conjecture 1]). Recently were obtained strong results 
in this direction (see [73], [7], [44]). However, generally the problem is still 
far from being solved even for the classical PA. We explain the main reason 
for this situation. The “strong asymptotics theory” is a far and deep gener¬ 
alization of the classical Bernstein-Szego asymptotic theory for the ordinary 
polynomials, orthogonal on the unit segment A := [—1,1] and related to the 
general case of non-Hermitian orthogonal polynomials (see [77], [10], [52, 
Section 1 and Section 6], [55], the paper by H. Widom [79], as well). 


'^Sometimes in the current paper we follow the Chudnovsky’s terminology [18]: under 
“Fade approximants” (and “Fade polynomials” respectively), we also mean the Hermite- 
Fade approximants, the multipoint Fade approximants, etc. 

®It is worth noting that seminal Stahl Theorem is much more general than the original 
conjecture of Nuttall. Stahl Theorem was proved under the condition that the set of 
singular points of the given multivalued function is of zero (logarithmic) capacity instead 
of being a finite set. 
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Given a germ® / of a multivalued analytic function / with a finite number 
of branch points, the seminal Stahl Theorem' gives a complete answer to 
the problem of limit zero-pole distribution of the classical PA to /. Stahl 
Theorem is of a quite general character. It not only admits multivalued 
functions with a finite number of branch points on the Riemann sphere, but 
also multivalued analytic functions with a singular set of zero (logarithmic) 
capacity. The keystone of Stahl theory [64]-[68] is the existence of a unique 
(up to a compact set of zero capacity) maximal domain for the given at the 
point z = oo multivalued function /. The so-called “maximal domain” of 
holomorphy of / is a domain D = D[f) 3 oo, such that the given germ / 
can be continued as holomorphic (analytic and single-valued) function from 
a neighborhood of the infinity point z = oo io D (the function / can be 
continued analytically along each path that belongs in to D). “Maximal” 
means that dD is of “minimal capacity” among all compact sets dG, such 
that G is a domain, oo € G and / G J^{G). Such “maximal” domain D is 
unique (up to a compact set of zero capacity). The compact set S = S{f) := 
dD is now called “Stahl’s compact set” or “Stahl’s 5-compact set”, and D is 
called “Stahl’s domain”. The crucial properties of 5 for the Stahl theory to 
be true are the following: the complement D = C \ 5 is a domain, 5 consists 
of a finite number of analytic arcs (in fact, the union of the closures of the 
critical trajectories of a quadratic differential), and, finally, 5 possesses the 
so-called property of “symmetry”®, that is. 


dgjj{z,oo) 

dn+ 


dg^{z,oo) 

dn~ 


zG S' 


0 


( 11 ) 


where gjoizjOo) is Green’s function for the domain D, with the logarithmic 
singularity at the point .2 = oo, 5® is the union of all open arcs of 5 (which 
closures constitute 5, that is 5\5® is a finite® set), and dn~^ and dn~ are the 
inner (with respect to D) normal derivatives of goiz, oo) at a point z G 5® 
from the opposite sides of 5®. 

Given a finite set S C C, < oo, let i2/(C\S) be the set of all functions 
analytic in the domain C \ S. Let \ S) := J2/(C \ S) \ Jif{C \ S), i.e. 

for a fixed S a function / is from the set \ S) if it is a multivalued 

analytic function in the domain C \ S, but not a holomorphic function in 
C \ S (/ is analytic, but not single-valued). 

Let S = Sy = {ai,..., Op}, = p < oo, be the set of all branch points 
of/, i.e. / G j:/(C \ S) \ J^(C \ S) =: .£/0(C\S). Glearly, if S/ = {-1,1}, 
then Stahl’s compact set (with respect to the infinity point) 5 = [—1,1]. 

The properties of the compact set 5 described above are crucial for the 
validness of Stahl Theorem. 

We recall some basic notations from the general PA-theory. 


®For example, a convergent power series at a point zq G C; commonly we set zq = oo 
or zo = 0. 

^Because of its general character and the various subjects of complex analysis involved 
into the proof of the theorem, it is sometimes considered as “Stahl Theory”. 

Q 

Compact sets of such type are usually referred to as “S'-compact sets” or “S'-curves”, 
see [42], [61], [59], [32]. 

®In the general Stahl Theorem the set S \ S'® is of zero capacity. 
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Let there be given a function / G \ S), which is holomorphic at 

the infinity point, / G J^(oo). To be more precise, suppose we are given a 
convergent power series, i.e. a germ, 


OO 

f = E 


k=0 


Ck 

zk+l ’ 


( 12 ) 


and the corresponding analytic function / G j 2 /^(C \ S), where < oo. 
For a positive Borel measure /r with a compact support supp(/i) d C, denote 
by V^{z) the logarithmic potential^® of fi, that is: 

V^^{z):=f log I ^ dniC). 

J supp /i I ^ ^1 

Given a number n G N, let^^ Prifl^Pn,! £ Cn{z), Pn,i ^ 0, be the Fade 
polynomials (at the infinity point) of the function /. We recall that 


Rn{z) 


{Pn,0 + Pn,lf){z) 



Z ^ OO. 


(13) 


The polynomials Pn,o, Pn,i are not unique, but the ratio Pn,o/Pn,i is uniquely 
determined by the relation (13). The rational function [n/n]f := —Pn,o/Pn,i 
is called the diagonal Fade approximant of order n (or the n-diagonal Fade 
approximant) of the function / (at the infinity point), and Rn is the remain¬ 
der function of the Fade approximant. 

Given an arbitrary polynomial Q G C( 2 ;), Q ^ 0, denote by 


X{Q):= 

C:Q(C)=0 


the associated zero counting measure of Q (as usual, 5^ denotes the Dirac 
measure concentrated at the point C £ C). 

One of the main results of Stahl theory [64]-[68] is 


Stahl Theorem (H. Stahl, see [67], [71]). Let f G =^(oo), / G 

< oo. Let D = D[f) he Stahl’s “maximal” domain for f, S = dD - 
Stahl’s eompact set, [n/n\f = —Pnfl/Pn,i - the n-diagonal Pade approximant 
to the function f. Then the following statements are valid: 

1) there exists LZD of Pade polynomials Pn,j, j = 0, 1, namely, 

1 * 

-X(-Pn,i)—A, as n->oo, j = 0,1, (14) 

where X = Xs is a unique probability equilibrium measure of the eompact 
set S; 

2) the n-diagonal Pade approximants converge in capacity to the func¬ 
tion f inside the domain D, 

[n/n]f{z) f{z), n —)■ oo, z € D-, (15) 


necessary, the potential should be spherically normalized, see [24], [14] and Re¬ 
mark 5. 

^^There are some ambiguities in this notations with (9), apologies for that. 






8 NIKOLAY R. IKONOMOV, RALITZA K. KOVACHEVA, AND SERGEY P. SUETIN 


3) the rate of the convergence in (15) is completely characterized by the 
equality 

\{f ~ ri ^ oo, z e D. (16) 


Remark 1. In item 1), Xg is the unique probability equilibrium measure of 
the compact set S, that is V^^{z) = 75 , z ^ S, where V^^{z) = — f log \z — 
Cl dXsiC) is the logarithmic potential of the measure A 5 , 75 is the Robin 
constant for S. The notation stands for convergence of measures 

in the weak-star topology. The notation in items 2) and 3) means 

convergence in capacity inside (on compact subsets of) the domain D. 

As it is basely known, gj:){z,oo) = 75 — ^'’^^(z). From here it follows 
immediately, that the 5-property (11) may be written in an equivalent form 
(cf. (25)), namely 


dV^ 

9n+ 


(^) 


^ cO 

^(z), Z€S. 


(17) 


Remark 2. Let Pn d ~ monic Fade polynomials. Then the 

statement of item 1 ) is equivalent to the following: 


\Pn,jiz)\^^"'^ n->oo, zGD, j = 0,1. 

Similarly, denote by Rn{z) the normalized remainder function, 


(18) 


Pniz) - ^n+l+in +■■■ ’ 

where in £ 1^, in = o{n) as n —>■ 00 . Then the statement of item 3) is 
equivalent to the following: 

\B*n{z)\^^'" ^ ^ n^oo, zeD. (19) 


Let Tjf = { — 1,1}. Clearly, in this 5 = [—1,1]. Hence, for the case of 
a general multivalued analytic function / G \ Sy), #S/ < 00 , the 

associated (with the function /) Stahl’s compact set 5 should be considered 
as a natural replacement of the unit segment A = [—1,1], which is the 
main object in the classical Bernstein-Szegd asymptotic theory. Recall that 
the classical Bernstein-Szegd formula of strong asymptotics for polynomials, 
orthogonal on the segment A with respect to a weight p defined on A and 
“smooth” enough, is presented in terms of Szegd function, associated with 
p, and the “mapping” function (p{z) = z + {z^ — 1)^^^, associated with A. 
Viewing Stahl Theorem, it is natural to conjecture, that in the formula of 
strong asymptotics of Fade polynomials Pnj Stahl’s compact set 5 should 
appear in a natural way, but not the segment A or the finite union of the 
real segments (see [10], [79], [52], [53, Conjecture 1], [7], [31, Conjecture 2 ]). 

The existence of Stahl’s 5-curve for an arbitrary multivalued analytic 
function with a finite set of singular points, on one hand, and Stahl The¬ 
orem of the LZD of Fade polynomials, on the other, should play a crucial 
role in the forthcoming development of the theory of the strong asymp¬ 
totics of Fade polynomials, associated with a multivalued analytic function, 
and the corresponding theory of strong convergence of Fade approximants 
(cf. [28], [7]). 
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It is worth noting that, in general, the property of a symmetry (in other 
words, S'-property of type (11), (25), see [61, formula (11)]) gives rise to a 
local analogue of complex conjugation operation and to the method of the 
construction of an associated two-sheeted or three-sheeted Riemann surface 
(see [71], [25], [61, § 6]). 

In 1987 A. A. Gonchar and E. A. Rakhmanov [24] extended the original 
Stahl’s conception of symmetry (the S-curve conception) to the conception 
of “symmetry in the presence of an external field”.The general LZD- 
theorem of orthogonal polynomials was proven mainly by using the new 
notions in [24]; i.e., for polynomials orthogonal with respect to a variable 
weight (depending on the degree of the polynomial, see also [22]). This 
new concept of the “presence of an external field” offered Gonchar and 
Rakhmanov the chance to solve the well known “1/9-problem” (see [78], 
also Wolfram MathWorld: One-Ninth Constant, it is also known as “Varga’s 
constant”, see Wolfram MathWorld: Varga’s Constant). Finally, Stahl’s 
and Gonchar-Rakhmanov’s results gave rise to the so-called GRS-method 
(Gonchar-Rakhmanov-Stahl-method), which solves the problem of LZD of 
generalized (non-Hermitian) quasi-orthogonal polynomials. The method it¬ 
self consists of three steps. 

First step: to state an appropriate theoretical potential extremal prob¬ 
lem associated with the initial problem of LZD of Fade polynomials^^, and 
then proving the existence of an “extremal” compact set in some family of 
“admissible” compact sets.^^ 

Second step: to prove that the extremal compact set possesses the special 
property of “symmetry”, i.e., the set is a “weighted” S-curve, associated 
with the problem under consideration. 

Third step: to establish that the LZD of the generalized orthogonal poly¬ 
nomials under question exists and coincides with the equilibrium probability 
measure (in the presence of the external field), which is concentrated on the 
weighted S'-curve. 

2.2. Strong asymptotics in classical Fade approximants theory. We 

give one of the possible forms of the strong asymptotics of Fade polynomials 
under the conditions of Stahl Theorem (cf. [31, § 2, Conjecture Ij). It 
is worth noting that in general (except for the case of genus zero) such a 
representation is not unique (see [19], [56], [39]). The reason is that there 
exist so-called spurious zeros of Fade polynomials, with behavior, as n —>■ oo, 
that can be described in many ways. 

Conjecture 1 (cf. [31], § 2, Conjecture 1). Let f G Jf(oo) and f G 
j^°(C\F) for some finite set S C C. Let D = D{f) be Stahl’s maximal do¬ 
main associated with f, S = S{f) = dD - the corresponding Stahl’s eompact 

^^See also the pioneering paper by Gonchar and Rakhmanov [21], 1981, where for the 
first time the notion of logarithmic theoretical potential problem with an external field 
was introduced in connection with the problem of LZD of type II HP polynomials for the 

Angelesco systems of Markov-type functions. 

1 

That is, Fade polynomials, multipoint Fade polynomials, Hermite-Pade polynomials, 
generalized Fade polynomials, etc., see [18]. 

^^Sometimes it is called a “stationary” compact set of the associated extremal theoretic 
potential problem, see [43], [61]. 
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set for f, 9^2 = 5 ^ 2 (/) “ ihe canonical^^ hyperelliptic two-sheeted Riemann 
surface associated^^ with the compact set S, z = = {z, ±) G 9^2 “ 

an arbitrary point on the two-sheeted 912? 'l'n(z) = /) - Nuttall’s psi- 

function associated ^ ' with f and 912 • Then, after a suitable normalization of 
the Fade polynomials Pnj{z) = Pn,j{z; f), j = 0, 1, and the remainder func¬ 
tion Rn, the following relations take place in capacity inside the domain D: 

Pn,j{z) + o(l)), n oo, 

Rn{z) =’’ (1 + o(l))> n-^oo, (20) 

where w'^ = H 2 g+ 2 {z), -f^ 2 g +2 £ C 23 + 2 N; is the equation that determines^^ 

the hyperelliptic Riemann surface 9 I 2 of genus g. 

Conjecture 1 is based partly on Nuttall’s result from 1986 [53] (see also 
[44] ) via Liouville-Steklov asymptotic method for the function / of the spe¬ 
cial form 

3 3 

f{z) = Yl{z-ajT\ ajGC\Z, /(oo) = 1 ( 21 ) 

j=i i=i 

(cf. [53, § 5. Asymptotic Conjecture]). 

Conjecture 1 was proven in some particular cases in [70], [73], [7], comp, 
also [52], [55], [ 6 ], [31]. 

In general, the problem of the strong asymptotics of Fade polynomials is 
still open and the state is as follows: there are only partial results about 
the subject, and all of them are based on the existence of Stahl’s com¬ 
pact set S and the concepts associated with it, namely, the two-sheeted 
hyperelliptic Riemann surface and the corresponding Abelian integrals on 
it. The problem does not depend on the method used to produce the asymp¬ 
totic formula: method of matrix Riemann-Hilbert boundary value problem, 
Liouville-Steklov method, method of the Nuttall’s singular integral equa¬ 
tion. Ultimately, all these methods are based on the existence of Stahl’s 
5-curve. 

3. Two-point Fade approximants (or T-fractions) convergence 

THEORY 

Recently V. I. Buslaev [14] (see also [17], [15]) applied the GRS-method in 
treating the multipoint^® FA of a multivalued analytic function with a finite 
set of branch points. Buslaev [14] proved an analogue of the classical Stahl 
Theorem, furthermore, he discovered some special features of the multipoint 
FA. In particular, for the case of two-point FA, say at the points z = 0 and 
z = 00 , he was the first, who found that in the “generic case” the weighted 

Stahl’s sense, see [71] and also [6], [31]. 

^®Such Riemann surface is uniquely determined by S. 

l^See [73], ]75], [6], [31]. 

^®Recall that all zeros of the polynomial H 2 g +2 are simple. 

^®To be more precise, m-point where m € Z>o is fixed. 
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S'-curve^®, associated with this problem, partitions the Riemann sphere in 
an “optimal way” into two domains Dq 3 0 and D^o 3 oo. At the same 
time, it was well known that the GRS-method does not work in such a dis¬ 
connected situation. In the case of multipoint (m-point) PA the situation is 
even more complicated, since it should be treated as a problem of an optimal 
partition of the Riemann sphere into m domains. Nevertheless, the suitable 
generalization and improvement of the original version of the GRS-method 
to this situation was given by Buslaev in [14]. This generalization provided 
Buslaev the opportunity to extend Stahl theory to multipoint PA for mul¬ 
tivalued analytic functions.In fact, Buslaev made a new step towards the 
development of the GRS-method in order to extend it to the disconnected 
case. Ultimately, it made the GRS-method much more powerful than it was 
before. 

3.1. Buslaev’s S-curve and Buslaev’s limit zero-pole distribution 
theorem in two-point Fade approximants theory. Let us discuss in 
short the case of multipoint PA and Buslaev Theorem. The original version 
of Buslaev Theorem is of general kind, since it deals with m-point PA. To 
be more precise, given a set of m distinct points G C and a set 

f = {/i,..., /m} of m germs of m analytic functions, such that fj G Jif{zj), 
j = 1,... ,m, we seek a rational function G Cn{z) in z of order ^ n, and 
such that the following relations^^ hold: 

{fj - Bn){z) = 0{{Z - Zj)^^), Z^Zj, (22) 

m 

where ^ = 2re -|- 1, rij G Z+, j = 1,..., m. Suppose that all functions fj 

are multivalued analytic functions with finite sets of branch points. Gener- 
ically, all functions fj are distinct, that is, not even one of fj is an analytic 
continuation of another, say fk, ^ / j- 

Suppose that each fj G s/^{C \ Sj) := .e/(C \ Sj) \ J^{C \ Sj), where 
Sj = T,f. and < oo, that is, each fj is a multivalued analytic function 
on the Riemann sphere, which is punctured at a finite set of branch points 
of fj. Suppose that in (22), preserving the conditions of Bnslaev Theorem, 
rij/n -3 2pj as n —>■ oo, where Yfjf=iPj — Ij Pj ^ 0- Buslaev Theorem 
states that there exists, in the generic case, an “optimal partition” of the 
Riemann sphere into m domains Dj 9 Zj , which are “centered” at the given 
interpolation points Zj, and such that the compact set P := C \ ur=i D, 
consists of a finite number of analytic curves and possesses the property of 
“symmetry” (see (25)). The compact set P itself is an 5-curve “weighted” 
in the external field, given by the logarithmic potential of the negative unit 

on 

It is natural to call such a curve a “Buslaev’s 5-compact set”. 

^^And even to extend the Stahl theory to more general situation, namely, when the 
m-point PA corresponds to a set of germs f = {/i,..., fm} = {(/i, Zi), ..., {fm, Zm)} of 
m multivalued analytic functions, each of which has a finite set of singular points on the 
Riemann sphere. These m germs are given at m distinct points zi,... ,Zm and it may 
happen that not even one of these germs can be obtained via the analytic continuation of 
another germ. 

^^The corresponding relation should be changed when Zj = oo for some j € {!,..., m}, 
see (27) below. 
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charge —p, concentrated at the set of the given points {zi ,..., Zm}- 

m 

u = (23) 

Furthermore, each function fj G Jif{zj) possesses a holomorphic (analytic 
and single-valued) continuation into Dj, fj G j = 1,... ,m. Finally, 

the m-point PA Bn converges in capacity, as n —)• oo, inside each domain 
Dj to the function fj, 

Bn{z) n —)■ OO, z G Dj. (24) 


There exists a limit zero-pole distribution of the multipoint PA Bn- The 
distribution in question coincides with the probability measure A = Ar, 
which is concentrated on F and is equilibrium in the external field V~^. 
This field is determined by the negative charge concentrated on the set 
{zi,..., Zm} of the interpolation points, each having the “weight” of —pj at 
each point Zj (see (23)). Also, the 5-property of the curve F and the rate of 
convergence in (24) may be completely characterized by (cf. (11) and (17)): 


d{V^ - V^) _ diV^ - V^) 

■(^) ~ E Z (^)> 


(9n+ 


dn 


zGT^ 


and (cf. (16)) 


(25) 


\fj{z) - Bn{z)\ 


1 /n cap 


-2 E PkQoA^^^k) 

e fc=i 


n 


oo, 


G D 




(26) 


where go,, {z, zt) is the Green’s function for the domain D^, with a singularity 
at the point z = Zk G Dk, k = 1,... ,m. 

In what follows, for the sake of simplicity, we restrict our attention to the 
particular case m = 2 of Buslaev Theorem. Thus, we will discuss in details 
only the case of two-point Fade approximant. 

Let zi = 0, Z 2 = oo and f = {/o,/oo} be the set of two multivalued 
analytic functions, such that /o G J^{0) and /oo G ,^(oo), and also /o, /oo £ 
\ S), where < oo. Thus, each of the functions /o and /oo is a 
multivalued analytic function on the Riemann sphere, punctured at a finite 
set of points, each point is a branch point of /o or of /oo or of both of them. 
In other words, /o and /oo are the two germs of the multivalued analytic 
functions, given at the points zi = 0 and Z 2 = oo, respectively. It is worth 
noting that they may be two germs of the same analytic function, taken at 
two different points, namely at z = 0 and at z = oo. 

The two-point^^ PA is defined as follows. Given n G N, let Pn, Qn £ Cn{z), 
Qn ^ 0 , be polynomials of degree ^ n, such that^^ the following relations 
hold 


Rniz) := {Qnf-Pn)iz) 


0{z^), z^O, 
0{l/z), z ^ oo. 


(27) 


QQ 

In the classical terminology, this is the n-th truncated fraction of the classical T- 
fraction, see [29], and also [13]. 

^^For a fixed n G N, we can also claim that the left side of (27) is as a —> 0 

and 0(1) as 2 —>■ oo, but this does not change the convergence theorem itself. 
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The pair of polynomials Pn and Qn is not unique, but the rational function 
Bn ■= PnIQn is Uniquely determined by (27), and is called the two-point 
diagonal PA to the set of germs of the functions f = {/o, /oo}- In the generic 
case, from (27) follows that 



z —)• 0, 


(28) 


z —)• oo. 


If it exists, the rational function Bn = Bn{z; f) G C„(z) is uniquely deter¬ 
mined by the relation (28). 

As for the classical Stahl’s case, the existence of an S'-curve, associated 
with the two-point PA and weighted in the external field V~^, v = {Sq + 
(5oo)/2, is the crucial element for Buslaev’s two-point convergence theorem. 
Such a weighted 5-curve P = rBus(/o) /oo) exists^^ (see [17]) and makes an 
“optimal” partition of the Riemann sphere into two domains Dq 3 0 and 
Doc 3 oo, such that C = T>o LI T U Doo, fo £ ^{Dq) and /oo G Jif{Doo)- 
The compact set T is a weighted 5-curve, i.e. T consists of a finite number 
of analytic arcs and possesses the following property of “symmetry”: 


where A = Ar is the probability measure concentrated on T and the equi¬ 
librium measure in the external field V~'^{z) = | log 1^1; furthermore, A is 
generated by the negative unit charge —u, that is. 



As before, T^ is the union of all open arcs of T (the closures of which consti¬ 


tute T) and and dn~ are the inner (with respect to Dq and Doo) normal 
derivatives at a point z G T*^ from the opposite sides of T^. Clearly, A is the 
balayage of the measure u from Dq U D^ onto T. It is worth noting that 
T itself is a union of the closures of the critical trajectories of a quadratic 
differential (see [17], [14], [33]). 

Here, for the sake of simplicity, we only consider the case of two-point PA, 
and we set zi = 0 and Z 2 = oo. In what follows, we also suppose that /o and 
/oo are the germs of the same multivalued analytic function /, and denote 
them by /o G J^(0) and /oo G Jif{oo). We suppose that the function / 
has a finite set of singular points in C. For example, / may be an algebraic 
function, i.e. a function given by an algebraic equation over the field C(2;), 
or a solution of a linear homogeneous differential equation with polynomial 
coefficients from the ring C(z) (see [35], [57], [18], [44]). 

Notice that the functions fQ{z) = (1 — ~ 1, z —>• 0, and /oo = 

(z^ — 1)“^/^ ^ 1/z, z —>■ oo, are the germs of the same analytic function /, 
given by the equation ( 2 ^ — l)w'^ = 1. But the functions /o(^) = (1 — 
and /oo = (^^ — 1)“^'^^ -|- 1 are not so. Thus, the latest case is the generic 
case, and hence Dq n D^o = 0 (see Fig. 5). 

Now we are ready to formulate the particular case of Buslaev Theorem 
for two-point PA (cf. Stahl Theorem). 


25 


'In general, there may exist some degenerated cases. 
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Buslaev Two-Point Theorem (V. I. Buslaev [14]). Let the function 
f G n J^{oo), f G \ S), < oo, and let the pair of germs 

/0)/oo in a common positior?^. Let C = Dq U F U Hoo the optimal 
partition of the Riemann sphere into two domains Dq 3 0 and Doo 3 oo, 
such that /o G /oo G Dq n Doo = 0, and F possesses the 

weighted S-property with respect to the external field V~'^, v = (do + (^oo)/2. 
Then for the n-diagonal two-point PA Bn of the set of the germs f = {/o, /oo} 
the following statements hold true: 

1) there exists a limit zero-pole distribution for Bn, namely, 

-x(-Pn),-x(Qn) ^ Ar, ra-^oo; (31) 

n n 

2) there is convergence in capacity as n ^ oo, namely, 

Bn{z)^fo{z), ZGDo, Bn{z) ^ foo{z), Z ^ Doo] (32) 

3) the rate of convergence in (32) as n ^ oo is completely characterized 
by the relations 

\foiz) - B„(z)|'/" ^ z e Do, 

\foo{z) - Bniz)\^^^ ^ e-^-oo(---), ^ e Doo. 

Remark 3. As in Stahl Theorem, item 1) is equivalent to the following 
relation as n —>■ oo for the monic two-point Fade polynomials: 

^ z€DoU Doo. (34) 

Remark 4. “Optimal” should be understood in connection with the S- 
property (29). In fact, this property means that the compact set F possesses 
a “stationary” or “equilibrium” property in the presence of the external field 
V~'^. This stationary property is of unstable type, see [59]. 

In Buslaev Two-Point Theorem, equality (29) is more complicated than 
equality (11) in Stahl Theorem. To be more precise, equality (29) should 
be understood as follows. In the generic case, the compact set F divides 
the complex plane into two domains Dq 3 0 and Doo 3 oo, such that all 
three sets 7 := ODq n dDoo, 70 := dDo \ dDoo and 700 := dDoo \ dDo are 
nonempty sets. When 2; G 7, dn~ is the normal derivative at the point z 
from the boundary to the inside of Dq and is the normal derivative at 
the point z from the boundary to the inside of Doo. 11-2 G 70, dn~ and 
dn'^ are the normal derivatives at the point z to the inside of Dq from the 
opposite sides of 70. We treat the case z G 700 in a similar way, but with 
respect to the domain Doo. 

It should be emphasized that the main problem in this direction is to 
prove the existence of a stationary compact set F, and to characterize it as a 
weighted 5-curve. The fact, that the equilibrium measure A is the balayage 
of the measure u from Dq U Doo onto F, is a trivial one. 

Remark 5. We regard the “optimal” partition as optimal with respect to 
the given field^^ V~'^, where v = ^(5o + ^oo)- In general, in the two-point 


Equivalently, we say that Buslaev’s S'-curve F divides the Riemann sphere into two 
domains. 

^'^Here the potential V~‘'(z) = | log| 2 | is spherically normalized. 
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case we have v = po^o+Poo^oo, wherepoiPoo > 0, po+Poo = 1- Thus, for the 
hxed function / G \ S) and the fixed points zq = {) and z^o = oo, the 

optimal partition and the compact set T depend on the pair po,Poo (see [14]). 

3.2. Strong asymptotics in two-point Fade approximants theory. 

It is worth noting that, in general, the GRS-method is far from being com¬ 
plete, because in the theory of LZD of HP polynomials some complicated 
theoretical potential problems arise in a natural way. At the moment these 
problems cannot be solved via Buslaev’s approach^® (see [60], [61, Sec 6.2]). 
The reason is that in the case of the m-point PA, the external field is gener¬ 
ated by the finite number of m positive pointed charges, concentrated at the 
m interpolation points (see also [12], [16]). In contrast to this fact, in the 
case of HP polynomials the external field is of more complicated structure 
(see for example [61]). 

In the case of multipoint PA (in part, in the case of two-point PA), the sit¬ 
uation is similar to the case of the classical PA. Namely, given Buslaev’s com¬ 
pact set (i.e. weighted 5-curve) Pbus, the associated canonical two-sheeted 
Riemann surface immediately appears, and it is equipped with the corre¬ 
sponding Abelian integrals, etc. Recently, A. V. Komlov and S. P. Suetin [33] 
derived a formula for the strong asymptotics of the two-point Pade poly¬ 
nomials for a special class of multivalued analytic functions, given by the 
representation 



(35) 


The result was obtained by use of the following reasoning: first, the authors 
proved that the corresponding two-point Pade polynomials of degree n and 
the corresponding remainder function are just the independent solutions of 
a linear differential equation of second order, which contains some poly¬ 
nomial accessory parameter of fixed degree, but depending on n. Second, 
they proved that Buslaev’s compact set, associated with the function (35), 
yields the Stokes lines^® for this differential equation. Finally, following 
Nuttall’s approach [53] (see also [44]), a suitable modification of the classi¬ 
cal Liouville-Steklov method (in accordance with the presence of accessory 
parameter) was used to derive formulae of strong asymptotics of two-point 
Pade polynomials®” and the remainder function. It follows from the for¬ 
mer and the latter, that the existence of Stahl’s 5-curve, in the case of the 
classical PA, and the existence of Buslaev’s weighted 5-curve, in the case of 
two-point PA, play the crucial role not only in the problem of limit zero-pole 
distribution of PA, but also in the problem of strong asymptotics of PA in 
both cases. Thus, for the moment, the problem of weak-star asymptotics 
is completely solved in both cases (for classical PA and for m-point PA it 
was completed by H. Stahl [64]-[68] and by V. I. Buslaev [14], [17], [15], 

QQ 

In Buslaev Theorem, a set f of germs of multivalued functions fj shonld satisfy a 
more restricted conditions, than in the original Stahl Theorem. Namely, f should have 
a finite set of singular points instead of a set of zero capacity, as it supposed in Stahl 
Theorem. 

on 

Since the accessory parameter depends on n, it would be better to say that Buslaev’s 
compact set attracts the Stokes lines of the differential equation, as n —>■ oo. 

®”lt is natural to call these polynomials two-point Jacobi polynomials, see (36). 
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respectively). These results should be regarded as successful applications 
of the GRS-method. Notice that the first result is connected with the pure 
logarithmic equilibrium problem, without any external field, and the second 
one is connected with logarithmic equilibrium problem under the presence 
of the external field. In general, the problem of strong asymptotics of the 
Fade polynomials is still open in both cases. However, for the moment, some 
promising results are obtained in both cases, and all of them are based on 
the existence of the associated 5-curve, and the associated weighted 5-curve 
as well. 

Notice that the two-point Fade polynomial Qn of the function /, given 
by (35), satishes the following non-Hermitian orthogonality relations, with 
respect to the variable weight function — ai)“(C “ ®2)~“/C"') 

where 7 is an arbitrary contour, that creates a path through the points 
ai, 02 , and separates the point z = 0 from the infinity point z = 00 . 

In the general theory of HF polynomials, the situation is quite different 
from above. In that direction, the hrst results of general character were 
obtained by A. A. Gonchar and E. A. Rakhmanov [21] in 1981 for the case 
of pure Angelesco systems of Markov-type functions. In 1984, J. Nuttall [52] 
stated some conjectures about HF polynomials, that were important and 
with an impact to the theory of HF polynomials. In 1986, E. M. Nikishin [48] 
discovered and investigated new systems of Markov-type functions, that are 
still of great interest for all experts in HF polynomials theory.^^ These 
systems are named after him, Nikishin systems. In 1988, H. Stahl [69] 
published a survey,®^ where he submitted a lot of conjectures about LZD 
of type I and type H HF polynomials, that were based, in part, on the 
paper by Gonchar and Rakhmanov [21]. As in Nuttall’s paper [52], the 
conjectures from [69] were about the real and the complex case. In 1997, 
A. A. Gonchar, E. A. Rakhmanov, and V. N. Sorokin [26] proved a general 
result about LZD of type H HF polynomials for the mixed, i.e. Angelesco 
and Nikishin, systems of Markov-type functions. In 2010, A. I. Aptekarev 
and V. G. Lysov [4] proved the most general at the moment result about LZD 
of type H HF polynomials for the mixed systems of Markov-type functions. 
To be more exact, they imposed the case when the support of one function 
of Markov type among the system under consideration is a proper part of the 
support of another Markov-type function among this system (for the proof 
of a partial case via another approach the reader is referred to [58]). All 
these general results [21], [26], [4] are dealing with systems of Markov-type 
functions with supports lying on the real line. But not even one of them 
may be applied to the situation under question in the current paper: when 
the supports of two Markov-type functions have nonempty intersection, but 
not even one of them is a proper part of the other. 


01 p 

One of the main reasons is that in the generic case the pair of functions /, / form a 
Nikishin system. 

^^The full manuscript of [69] was available much earlier than 1988. It is much bigger 
than the length of the published paper, and is now available in electronic form. 
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In the fundamental paper by Nuttall [52] were proposed several conjec¬ 
tures of a general type about asymptotics of HP polynomials. All of them 
are concerned with both type I and type II HP polynomials, and only with 
the strong asymptotics of HP polynomials. Neither the existence of the 
associated S-curve, nor the weak-star asymptotics of HP polynomials were 
discussed in [52]. This situation is typical for all results of HP polynomials. 
In fact, the only results of general character on this subject were obtained 
for the “real case” situation. In other words, all rigorous and general results 
in this direction were proved for the systems of Markov-type functions, and 
under such the conditions that LZD of HP is completely described by the 
associated “matrix of interaction” and extremal theoretical potential prob¬ 
lem. Depending on the system of Markov-type functions, this matrix may 
be of Angelesco type, or of Nikishin type, or some kind of “mixed” type. 
In each case, the problem of LZD of HP polynomials is solved in terms of 
the corresponding equilibrium measures, concentrated on finite number of 
segments of the real line. The most general description of such “real case” 
was done by Gonchar and Rakhmanov [23] (see also the above cited pa¬ 
pers [21], [48], [26], [4]). It is worth noting that some special properties of 
the HP polynomials in the real case were proved by G. Lopez Lagomasino 
and coauthors [38] , [37] , [40] , namely, the property of normality, the conver¬ 
gence property (but without any characterization of the rate of convergence), 
the interpolation property, etc. 

Thus, in all cases discussed above, the answer to the problem of LZD of 
HP polynomials was given in terms in equilibrium measures concentrated 
on the segments of the real line.^^ Presently, only a few rigorous results 
are known of the asymptotics of HP polynomials for the so-called “complex 
case” situation. In other words, if the solution of the problem of LZD 
of HP polynomials is given in terms of equilibrium measures concentrated 
on the associated S'-curves, then the measures are located somewhere in 
the complex plane, but not on the real line. As usual, in such situations, 
we should first prove the existence of an associated 5-curve^^ (cf. Stahl’s 
and Buslaev’s theorems). Also, only a few rigorous results of LZD of HP 
polynomials are known for the case when m = 3, that is, for the collection of 
three functions [1, /i, / 2 ] instead of the collection of two functions [1, /], as in 
the classical PA case. Notice that, for the moment, there are two different 
approaches to the problem under question. The first one is based on the 
cubic equation, and the other is based on the concept of Nuttall’s condenser. 
Here, we do not discuss the details, but instead refer the reader to the 
papers [54], [2], [8], [9], [3] and the references therein. For the cubic equation 
and Nuttall’s condenser, we refer to [30], [1], [2], [5], [46] and [60], [61], 
[76], [31], respectively. 

We can draw the following conclusion from the results of the papers re¬ 
viewed above. So far, there does not exist a general approach to the problem 

OQ 

Because of this, we refer to such type of case as “real case”. 

contrast to Stahl’s and Buslaev’s theorems, in the case of HP for the collection of 
three functions [1, /i, A], the associated S'-compact set should be considered as consisting 
of two proper subsets, that play different roles in the associated theoretical potential 
problem. The example of such S'-compact set is given by the “Nuttall’s condenser”, 
see [60], [61], [31]. 
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of LZD of HP polynomials, and there is no connection to the powerful GRS- 
method, and, moreover, there does not even exist a suitable conjecture on 
this subject, similar to Stahl’s and Buslaev’s results. In our knowledge, all 
rigorously proved results are of partial character, and there does not exist 
any theorem on LZD of HP polynomials, with which we can explain from 
a theoretical viewpoint the results of the numerical experiments, submitted 
in the current paper. The numerical results, presented herein, can be de¬ 
scribed in simple terms, but it is difficult to explain them from a theoretical 
viewpoint. Ultimately, these results can be regarded as a challenge to all 
experts on HP theory. 


4. Numerical results 

For completeness and for the reader’s convenience, at first we present 
some numerical results, concerned with Stahl’s and Buslaev’s theorems. 

4.1. Some numerical examples for classical Fade approximants. 

4.1.1. The case of a function with three branch points, Chebotarev-Stahl’s 
S-curve and Stahl’s limit zero-pole distribution theorem. Let 

fiz) = ^/[{z - ai){z - a 2 ){z - 03 ))^''^, (37) 

where oi = —1.2 -|- 0.8i, 02 = 0.9 -|- 1.5i, 03 = 0.5 — 1.2i. On Figure 1, 
the zeros (blue points) and the poles (red points) of the PA [130/130]/ of 
/ at infinity are plotted. On Fig. 3 and Fig. 2 the zeros (blue points) and 
the poles (red points) of PA [130/130]/ are plotted separately. Clearly, 
numerical zero-pole distribution is in good agreement with the statements 
of Stahl Theorem. But there is a spurious zero-pole pair, which does not 
correspond to any singularity of the given function (37). The behavior of 
this pair as n —t- 00 is not governed by Stahl Theorem, since this theorem 
is dealing with a weak-star limit zero-pole distribution of PA. This kind of 
behavior was discovered by Nuttall [53] in 1986 via a suitable modification 
of Liouville-Steklov method (see also [44]). However, Nuttall’s result was 
found only after Stahl Theorem was proved. To be more precise, Nuttall 
used the existence of Stahl’s S'-curve for the function (37). 

Namely, Nuttall first proved that Pade polynomials of order n and the 
remainder function for the function (37) both solve a linear homogeneous 
differential equation of second order, with polynomial coefficients of fixed 
degrees. These polynomials have some accessory parameters depending on n. 
Subsequently, Nuttall proved that Stahl’s 5-curve is the limit, as re —)• 00 , of 
the Stokes lines for these differential equations, that depend on re. Finally, 
in [53] Nuttall proved that there may exist only one spurious zero-pole pair 
of PA to the function (37), since the two-sheeted Stahl’s Riemann surface 
y\. 2 , associated with the function (37), is of genus g = 1. Nuttall proved that 
if the zero and the pole of a spurious zero-pole pair are close to one another, 
then they actually cancel each other, as re —)> 00 . But if, for each re G N, 
they do not coincide with each other, then they are everywhere dense on the 
Riemann sphere, as re —)• 00 . Describing the behavior of the spurious zero- 
pole pairs might require solving the so-called “Jacobi inversion problem” 
(see [63]), that is an equation with an Abelian integral of first kind on the 
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left and some expression that is linear in n on the right. Thus, as n —>■ oo, 
such pairs form some type of “winding of the torus”, which is everywhere 
dense on that torus, that is, dense on the Stahl’s two-sheeted Riemann 
surface ^2 of genus 5 = 1. The fact that the zero and the pole in some sense 
cancel each other in such a zero-pole pair, as n 00 , is in full agreement 
with the pure numerical results that were obtained by M. Froissart [20], see 
also [11, Chapter 2, §2.2] and [74]. 

4.1.2. The case of a function with several branch points, Chebotarev-Stahl’s 
S-curve and Stahl limit zero-pole distribution theorem. Let 



(38) 


where oi = 4.3 i, 02 = 2 -\- 0.5i, 03 = 2 -|- 2i, 04 = 1 — 3z, 05 = 4 -|- 2i, 
oq = 3-I-5L On Fig. 4 are plotted zeros (blue points) and poles (red points) of 
the PA [103/103]j, taken at the infinity point z = 00 to the function / given 
by (38). Clearly, numerical zero-pole distribution is in good agreement with 
the statements of Stahl Theorem. But there are spurious zero-pole pairs, 
which do not correspond to singularities of the given function (38). The be¬ 
havior of this pair, as n —)• 00, is not governed by Stahl Theorem, since it is 
about weak-star limit zero-pole distribution of PA. This behavior was discov¬ 
ered by A. Martinez-Finkelshtein, E. A. Rakhmanov, and S. P. Suetin [44] 
in 2012 via a suitable modification of Liouville-Steklov method (cf. [53]). 
This result was obtained on the base of Stahl Theorem. More precisely: the 
authors of [44] used the basic fact about the existence of Stahl’s 5-curve 
for the function (38). Namely, first they proved that Pade polynomials of 
order n, as well as the remainder function for the function (38) solve a linear 
homogeneous differential equation of second order with polynomial coeffi¬ 
cients of fixed degrees. These polynomials have some accessory parameters 
depending on n. Second, they showed that Stahl’s 5-curve is the limit, as 
n —)• 00, of the Stokes lines for these differential equations, depending on n. 
Ultimately, the authors concluded that no more than four spurious zero-pole 
pairs of PA of the function (38) may exist. The reason is that the Stahl’s 
two-sheeted Riemann surface 942 , which corresponds to the function (38), is 
of genus 5 = 4. It was proved that the zero and the pole in each spurious 
zero-pole pair are close to each other and they cancel each other, as n —>• 00. 
However, in the “generic case”, if they do not coincide with each other for 
each n G N, then they are everywhere dense on the Riemann sphere, as 
n —?• 00. Describing the behavior of the spurious zero-pole pairs might re¬ 
quire solving a system of equations with some Abelian integrals of first kind 
on the left side and some expressions, which are linear in n, on the right 
side. The fact that the zero and pole in some sense cancel each other in 
such a zero-pole pair as n —)• 00 is in full agreement with the pure numerical 
results obtained by M. Froissart [20], see also [11, Chapter 2, §2.2] and [74]. 

4.2. Some numerical examples for two-point Pade approximants. 

Let 



( 39 ) 
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where oi = 0.9 — l.lz, 02 = 0.1 + 0.2i. Suppose that two “different” germs 
are taken, /o and foo of the function (39) at the point z = 0 and at the 
infinity point z = 00 , respectively. Here, “different” means that to obtain 
the germ f^o £ .;^(oo) by the analytic continuation of the germ /o G =^( 0 ), 
we should go along a path that encircles exactly one time one of the branch 
points (oi or 02 ). 

On Fig. 6 are plotted the zeros and poles of two-point PA [199/199]/|3 
to the function (39). On Fig. 7 and Fig. 8 are plotted separately the zeros 
(blue points) and poles (red points) of the two-point PA [199/199]/g for 
the function (39). This numerical zero-pole distribution of two-point PA 
[199/199]/p is in full agreement with Buslaev’s limit zero-pole distribution 
theorem [14]. Similarly to the function (37) and classical PA, the two-sheeted 
Riemann surface, associated with the function (39) and the two-point PA, 
is of genus g = 1. Therefore, similarly to the case of classical PA, there is 
a single zero-pole pair of spurious character. The behavior of this pair, as 
n —?• 00 , cannot be described by Buslaev Theorem, since the theorem is only 
about weak-star convergence of two-point PA. 

The description of this behavior for the special function (39) and under 
the condition of a “generic case”, imposed on the branch points ai and 
02 , was done by Komlov and Suetin in [33]. Under this condition, in [33] 
was derived a formula of the strong asymptotics for two-point PA. The 
method producing this formula is similar to Nuttall’s method for classical 
PA, see [53] and [44]. In [33], it was proved first that for each n G N the two- 
point Pade polynomials and the corresponding remainder function solve the 
homogeneous linear differential equation of second order with polynomial 
coefficients of fixed degrees, but depending on n. After that, it was proved 
that Buslaev’s S'-curve forms the limit, as n —>• 00 , of the Stokes lines for 
these differential equations. Finally, a formula of the strong asymptotics for 
the two-point PA was found via the asymptotic Liouville-Steklov method. 

4.3. Some numerical examples for Hermite-Pade polynomials. Now 

we are ready to make some short description of the obtained numerical 
results for the HP of the collection of three functions [1, /i, / 2 ]. 

4.3.1. Parameter a = —0.1.' logarithmic functions, Fig. 9-Fig. 10, square 
root functions. Fig. 11-Fig. 12. At first we set a = —0.1 (and keep n = 200). 
Clearly, that value of a is out of the range of (0,1). But in that case all three 
pairs /i,/ 2 , 51,52 and /ii,/i 2 form Angelesco systems, since Ei,E 2 C M and 
EiD E 2 = 0. On Fig. 9, Fig. 10, Fig. 11, and Fig. 12 are plotted the zeros 
of HP polynomials Q 200 ,j and P 200 J, j = 0,1, 2. The numerical distribution 
of the zeros of the HP polynomials are similar in both cases, and it does not 
depend on the different type of branching of the functions fi , /2 and 51,52 at 
the points 2 ; = ±1, ±a. The reason is as follows. All these numerical results 
are in a full agreement with the general theory of limit zero distribution 
(LZD) for HP polynomials in the Angelesco case (see [21], [26], [4]; in these 
papers the case of type H HP polynomials was considered, but it is well 
known that in the Angelesco case there is a direct connection between LZD of 
type I and type H HP polynomials). Since the size of the supports of Markov- 
type functions fi and /2 (and 51 , 52 as well) are equal, the phenomena of 
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the so-called “pushing of the charge”^^ is absent in the case when a = —0.1. 
It is also in full agreement with the general theory (see [30], [1], [21]). From 
Fig. 9 and Fig. 10 follows, that the zeros of HP polynomials <5200,1 and 
< 5200,2 (and F^oo,!) ^ 200,2 as well) are located on the segments Ei and £^ 2 , 
respectively. But the zeros of < 5200,0 (and £ 200,0 respectively) are located 
on the imaginary axis. Thus, in some sense, the zeros of < 5200,0 “intend to 
separate” the zeros of < 5200,1 and < 5200,2 from each other. Again, this result is 
in good agreement with the well known LZD of type I HP polynomials in the 
Angelesco case (see, e.g. [54]). To hnish the analysis of the case a = —0.1, 
we emphasize that in this case all numerical results are in good agreement 
with the general theory of HP polynomials. Thus it is reasonable to consider 
all numerical results for a = 0.2; 0.4; 0.625; 0.73; 0.8 as to be trustable. 

4.3.2. Parameter a = 0.2; logarithmic functions, Fig. 13-Fig. 16. Now let 
a = 0.2. For this value of the parameter a, the intersection of the two 
segments £1 = [—1, a] and £2 = [—a, 1] equals £1 n £2 = A = [—a, a] 0, 
that is, A is small, but nonempty real segment. 

In all three cases, Case I, Case 2, and Case 3, the numerical zero distri¬ 
bution (NZD) of HP Q 200 ,j{z] fi, 12 ), P 200 ,jiz;gi,g 2 ), and U 2 ooj{z] hi, 62 ) 
for j = 1,2 are similar to each other. Thus, here we only discuss the NZD 
for HP < 5200,1 and < 5200,2 (these zeros are plotted by red and black points 
respectively, see Fig. 15 and Fig. 16). Notice that zeros of < 5200,1 and < 5200,2 
are symmetric to each other, with respect to the imaginary axis. Thus, it is 
enough to only consider the NZD for HP < 5200,2 (these zeros are plotted in 
black; see Fig. 16). 

For this value of the parameter a, the essential part of zeros of HP <5200,2 
are located not on the real line, but somewhere in the complex plane. Evi¬ 
dently, this set is symmetric with respect to the real axis, since we consider 
Markov-type functions. Due to the general conjectures (see [52], [69]), that 
are based on the seminal paper by Gonchar and Rakhmanov [21] about 
LZD of pure Angelesco system of Markov-type functions, the LZD for HP 
<5n,2, as n —>• 00, should be described via an equilibrium measure, say A 2 , 
that is associated with some special max-min extremal theoretical poten¬ 
tial problem. From Fig. 16 (see also Fig. 15) it follows immediately, that 
in the case a = 0.2 the support of the measure A 2 should be a discon¬ 
nected set. The same is true for the equilibrium measure Ai, that is as¬ 
sociated with LZD for HP <5n,i. The union S'i(/i,/ 2 ) U «S' 2 (/i,/ 2 ) of two 
compact sets S'i(/i,/ 2 ) = suppAi and 5'2(/i,/2) = suppA 2 form some type 
of lenses. From the general approach, based on the paper by Gonchar and 
Rakhmanov [24] about ‘T/9-conjecture”, it follows, that the compact sets 
Siifi, f 2 ) and 5'2 (/i,/ 2) should be a weighted 5-curve (see also [59]). No¬ 
tice, that in both cases the open sets C \ 5i(/i,/2) and C \ 52(/i,/2) are 
also domains (see Fig. 15 and Fig. 16 respectively). 

It follows from Fig. 14, that the LZD for HP < 5 n,o should be quite different 
from the LZD for HP < 5 n,i and < 5 n, 2 - The reason is as follows. This LZD 
should be described by an extremal max-min theoretical potential problem of 
different type than before. In part, the support 5o(/i, / 2 ) of the equilibrium 

^^The phenomena was discovered by Gonchar and Rakhmanov in [21], 1981. 
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measure Aq, associated with this extremal theoretical potential problem, 
should be a continuum, i.e. a connected compact set. But now, the open 
set C \ 5o(/i,/ 2) is not a domain, since it consists of three domains. The 
compact set S'o(/i, / 2 ) should also be a weighted 5-curve, but of some other 
nature than the compact sets 5i(/i,/2) and 52(/i,/2). In particular, there 
should be three Chebotarev’s points on the compact set 5o(/i,/2), i.e. the 
points of zero density of the equilibrium measure Aq. One of these points 
is located on the upper half plane, the other point is located on lower half 
plane, and the third Chebotarev’s point coincides with the infinity point. 
Notice, that similarly to the case when a = —0.1, now the 5-curve 5o(/i, f 2 ) 
separates the 5-curves 5i(/i, / 2 ) and 52(/i, / 2 ) from each other. It might be 
conjectured that the pair 5(/i,/ 2 ), 52(/i,/ 2 ) forms some kind of weighted 
Nuttall’s condenser (see [61]). 

4.3.3. Parameter a = O.2.- square root functions, Fig. 17-Fig. 20, and cubic 
root functions, Fig. 21-Fig. 2f. The NZD of HP P 200 ,i and 5*200,2 for the 
square root functions, given by (5), (6), and NZD of HP polynomials 1 / 200,1 
and 1 / 200,2 for the cubic root functions, given by (7), (8), are similar to the 
NZD of HP Q 200 ,i and Q 200 , 2 , for the pair of logarithmic functions given 
by (3), (4), respectively (see Fig. 19, Fig. 20, Fig. 23, and Fig. 24). It is also 
valid for the HP t/ 200 , 0 ) be. this NZD is similar to the NZD of HP Q 200,0 for 
the given logarithmic functions. Therefore, we do not discuss specially the 
NZD of HP 5*200,1, 5*200,2, t/ 200 , 1 , 1 / 200 , 2 , and 1 / 200 , 0 - Notice, that there is a 
pair of spurious zeros of the HP 1 / 200,0 (see Fig. 22). From that numerical 
fact follows, that the genus of the associated three-sheeted Riemann surface 
933 should be equal to 2 (cf. [54], [2]). 

4.3.4. Parameter a = 0.2; square root functions, NZD of 5*2oo,o, Fig. 18. 
The NZD of HP polynomials 5*2oo,o is quite different from NZD of HP Q 200,0 
(see Fig. 14) and 1 / 200,0 (see Fig. 22). In fact, from the numerical results 
follows, that the associated 5-curve So{gi,g 2 ) should consist of two seg¬ 
ments, and the open set C \ 5o(5i, 52 ) should be a domain. In that case, the 
compact set 50 ( 51 ,( 72 ) separates distinctly the other 5-curves 51 ( 51 , 52 ) and 
< 52 ( 51 , 52 ) from each other. 

4.3.5. Parameter a = 0.4; logarithmic functions. Fig. 25-Fig. 28, square 
root functions. Fig. 29-Fig. 32, cubic root functions. Fig. 33-Fig. 36; pa¬ 
rameter a = 0.625; logarithmic functions. Fig. 37-Fig. 40 , square root func¬ 
tions, Fig. 41 -Fig. 44} cubic root functions. Fig. 45-Fig. 48 . 

For these two values of the parameter a all the NZD of HP <5200,^, 5*200,i, 
1 / 200 ,^, j = 0,1,2, are similar to the case when a = 0.2. The observed 
difference does not change the principal structure of the associated 5-curves 
and is the following. The lenses, that appeared for parameter a = 0.2, 
become larger and larger when a changes from 0.2 to the values a = 0.4 and 
a = 0.625. The positions of all four vertices of the lenses are not fixed, but 
depend on the value of the parameter a. When the parameter a increases 
from a = 0.2 to a = 0.4 and after that to o = 0.625, the two real vertices 
move from the inside of the segment [—1,1] towards to the end points ±1. In 
addition, the two pure imaginary vertices move along simultaneously from 
the imaginary axis to the infinity point, where they meet each other, under 
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some critical value of the parameter a* G (0.625,0.73) (see the next Fig. 49- 
Fig. 72; the statement on the existence of the critical value a* is based only 
on numerical results and should be considered as a conjecture). As usually, 
the NZD of HP i- 200,0 corresponds to the whole segment [—a, a]. Hence, the 
blue segment on Fig. 18, Fig. 30, and Fig. 42 becomes wider and wider as 
the parameter a increases form a = 0.2 to a = 0.625. 

4.3.6. Parameter a = 0.73." logarithmie functions, Fig. 49-Fig. 52, square 
root functions, Fig. 53-Fig. 56, cubic root functions, Fig. 57-Fig. 60; param¬ 
eter a = 0.8." logarithmic functions. Fig. 61-Fig. 64 , square root functions. 
Fig. 65-Fig. 68, cubic root functions. Fig. 69-Fig. 72. After a> a*, where 
a* is the critical value of the parameter a described above, the NZD of HP 
Q 200 JJ -P 200 ,j) b^ 200 j) j = 0,1, 2, dramatically changes. Namely, the comple¬ 
ments C \ S'! and C \ 52 of the 5-curves 5i and S 2 are now disconnected open 
sets, instead of domains, as they were when a < a*. Plotted on a single pic¬ 
ture, all three numerical sets 5i,52,53, that consist of red, black, and blue 
points, respectively, form in these three cases three quite different structures. 
For instance, when the parameter a equals 0.8, in Case 1 the associated 
5-curve 5o(/i,/2) becomes a disconnected compact set, namely, 5o(/i,/2) 
consists of two continua (see Fig. 62). In Case 2, we have So{gi,g 2 ) = [—a, a] 
(see Fig. 66), and finally, in Case 3 the 5-curve 5o(/ii, 62 ) is a continua that 
contains the infinity point (see Fig. 70). In the last case, 5o(/ii,/i2) con¬ 
tains two Chebotarev’s points that are the points of zero density for the 
corresponding equilibrium measure Aq with suppAo = 5o(/ii,/i2). 
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Figure 1. Zeros and poles of the diagonal Fade approximant 
[130/130]/ of the function /(z) = l/((z— (—1.2-|-0.8 j)){z — (0.9-1- 
1.5i))(z — (0.5 — 1.2i)))^/^, distributed accordingly to the electro- 
statical model by E. A. Rakhmanov [59]. There is a Froissart 
doublet (spurious zero-pole pair) when n = 130 (see also Fig. 2 
and Fig. 3). Since the genus of the Riemann surface is 1, there 
might be at most one Froissart doublet. In full compliance with 
the Rakhmanov model [59], the Froissart doublet “attracts” the 
Stahl S'-compact iSiao- 
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Figure 2. The poles of the Fade approximant [130/130]/ approx¬ 
imate a Chebotarev point z;i 3 o for the S'-compact Siao (see [59]). 
The Chebotarev point is at (0.029,0.466). When n —>• oo we have 
that Vn ^ V is a classical Chebotarev point. There is one spurious 
pole of the Fade approximant [130/130]/, it is accompanied by a 
spurious zero of the Fade approximant [130/130]/ (see Fig. 3). The 
spurious zero-pole is at (0.469,0.633). 
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Figure 3. The Chebotarev point should not be approximated by 
zeros of the Fade approximant [130/130]/ of the function (37). 
Evidently, the Cheboratev point does not exist on the picture. 
There is one spurious zero of the Fade approximant [130/130]/, 
it is accompanied by a spurious pole of the Fade approximant 
[130/130]/ (see Fig. 2). 
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Figure 4. Zeros and poles of the diagonal Fade approximant 
[103/103]/ of the function f{z) = l/((z + (4.3 + 1.0i))(2r — (2.0 + 
0.5i))(z + (2.0 + 2.0i))iz + (1.0 - 3.0i))(z - (4.0 + 2.0i)){z - (3.0 + 
5.0i)))^/®. These zeros and poles are distributed in a plane, un¬ 
der fixed n = 103, accordingly to the electrostatical model by 
Rakhmanov [59]. Since the genus of the Riemann surface is 4, 
for each n there might be no more than 4 Froissart doublets. Here 
are observed 3 Froissart doublets. In full compliance with the 
Rakhmanov model, the Froissart doublets “attract” the Stahl S- 
compact S'io 3 - In general, the zeros and poles of the diagonal 
Fade approximants [n/n]f are distributed as n —>■ oo accordingly 
to Stahl Theorem [71]. 
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Figure 5. Numerical zeros and poles distribution of two-point 
Fade approximants [120/120] to the set of functions f = {/o, /oo}, 
where fo = ((1 - 2z){2 - - l){z - 2))-'^/'^ + 1. 

The germs /g and f^o result in two different multivalued analytic 
functions. Thus, this is a generic case and by Buslaev Theorem 
the associated S-curve partitions the Riemann sphere into two do¬ 


mains. 
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Figure 6. Numerical zeros and poles distribution of two-point 
Fade approximants [199/199] to the function f{z) = {{z — ai)/{z — 
02 ))^/"^, where Oi = 0.9 — l.li and 02 = 0.1 -I- 0.2i. Here are 
selected two “different branches” of the function /, namely, /o = 
((z —ai)/( 2 ; — 02 ))^/"* and foo = —{{z — ai)/{z — a 2 ))^^^. Almost all 
zeros (blue points) and poles (red points) approximate numerically 
Buslaev’s compact set. But there is a spurious zero-pole pair that 
moves as n —> 00 under the order of equation from [33]. 
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Figure 7. Numerical zeros distribution of two-point Fade approx- 
imants [199/199] to the function f{z) = {{z — ai)/{z — 02 ))^^"^, 
where ai = 0.9 — 1.1* and 02 = 0.1 + 0.2*. Here are se¬ 
lected two “different branches” of the function /, namely, /o = 
{{z - ai)/{z - 02 ))^/"* and foe = -{{z - ai)/{z - 02 ))^^^. Almost 
all zeros (blue points) approximate numerically Buslaev’s compact 
set. But there is a spurious zero-pole pair that moves as n —)■ 00 
under the order of equation from [33]. 
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Figure 8. Numerical poles distribution of two-point Fade approx- 
imants [199/199] to the function f{z) = {{z — ai)/{z — 02 ))^^"^, 
where ai = 0.9 — 1.1* and 02 = 0.1 + 0.2*. Here are se¬ 
lected two “different branches” of the function /, namely, /o = 
{{z - ai)/{z - 02 ))^/"* and foo = -{{z - ai)/{z - 02 ))^^^. Almost 
all poles (red points) approximate numerically Buslaev’s compact 
set. But there is a spurious zero-pole pair that moves as n —)■ 00 
under the order of equation from [33]. 
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Figure 9. Numerical distribution of zeros of type I Hermite- 
Pade polynomials < 5200,0 (blue points), <5200,1 (red points), <5200,2 
(black points), for the collection of functions [l,/i,/2], where 
/i = log((0.1 -b l/z)/(l + l/z)), /2 = log((0.1 - l/z)/{l - 1 /z)). 
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Figure 10. Numerical distribution of zeros of type I Hermite- 
Pade polynomials < 5200,0 (blue points) for the collection of func¬ 
tions [l,/i,/ 2 ], where /i = log((0.1 + l/ 2 :)/(l -b 1/z)), /2 = 
log((0.1-l/z)/(l-l/z)). 
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Figure 11. Numerical distribution of zeros of type I Hermite- 
Pade polynomials ^200,0 (blue points), ^200,1 (red points), ^200,2 
(black points), for the collection of functions [1,31,52], where 31 = 
((0.1 -f l/z)/(l + \lz)Yl\ 32 = ((0.1 - l/z)/(l - l/z))i/2. 
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Figure 12. Numerical distribution of zeros of type I Hermite- 
Pade polynomials ^200,0 (blue points) for the collection of functions 
[1,31,32], where 31 = ((0.1 -b llz)/{l + l/2))^^^ 92 = ((0.1 - 
l/z)l{l-l/z)f/\ 
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Figure 13. Numerical distribution of zeros of type I Hermite- 
Pade polynomials < 5200,0 (blue points), <5200,1 (red points), <5200,2 
(black points), for the collection of functions [1, /i, / 2 ], where /i = 
log((0.2 - l/ 2 :)/(l + l/z)), /2 = log((0.2 -b l/2;)/(l - l/z)). 
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Figure 14. Numerical distribution of zeros of type I Hermite- 
Pade polynomials <5200,0 (blue points) for the collection of func¬ 
tions [l,/i,/ 2 ], where /i = log((0.2 - l/z)/{I + l/z)), /2 = 
log((0.2 + l/z)/(l-l/z)). 
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Figure 15. Numerical distribution of zeros of type I Hermite- 
Pade polynomials Q 200 ,i (red points) for the collection of functions 
[l,/i,/ 2 ], where /i = iog((0.2 - l/z)/{l + l/z)), /2 = log((0.2 + 
l/z)/{l-l/z)). 



Figure 16. Numerical distribution of zeros of type I Hermite- 
Pade polynomials ( 5200,2 (black points) for the collection of func¬ 
tions [l,/i,/ 2 ], where fi = log((0.2 - l/z)/(l -b 1/z)), /2 = 
log((0.2 + l/z)/(l-l/z)). 
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Figure 17. Numerical distribution of zeros of type I Hermite- 
Pade polynomials F 200,0 (blue points), F 200 ,i (red points), ^ 200,2 
(black points), for the collection of functions [ 1 , 31 , 52 ], where 31 = 
((0.2 - l/z)/{l + l/z)Yl\ 32 = ((0.2 + l/z)/(l - l/z))i/2. 
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Figure 18. Numerical distribution of zeros of type I Hermite- 
Pade polynomials P 2 oo,o (blue points) for the collection of functions 
[ 1 , 51 , 52 ], where 31 = ((0.2 - llz)/{l + 1 / 2 ))^^^, 32 = ((0.2 -t 
l/z)l{l-l/z)f/\ 
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Figure 19. Numerical distribution of zeros of type I Hermite- 
Pade polynomials P 200.1 (red points) for the collection of functions 
[ 1 , 51 , 52 ], where gi = ((0.2 - l/z)/(l + l/z))^/^ 52 = ((0.2 + 



Figure 20. Numerical distribution of zeros of type I Hermite- 
Pade polynomials P200,2 (black points) for the collection of func¬ 
tions [1,51,52], where 51 = ((0.2 - l/z)/{l + llz))^/"^, 52 = 

((0.2 + l/z)/(l-l/z))V2. 
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Figure 21. Numerical distribution of zeros of type I Hermite- 
Pade polynomials U 2 oo,o (blue points), 1 / 200,1 (red points), 1 / 200,2 
(black points), for the collection of functions [1, hi, I 12 ]) where hi = 
((0.2 - l/z)/{l + h 2 = ((0.2 + l/z)/(l - l/z))i/ 3 . 
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Figure 22. Numerical distribution of zeros of type I Hermite- 
Pade polynomials 1 / 200,0 (blue points) for the collection of functions 
[l,hi,h 2 ], where hi = ((0.2 - l/z)/{l + 1/z))^/^, li 2 = ((0.2 -|- 
l/z)/{l-l/z)Y/\ 
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Figure 23. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,1 (red points) for the collection of functions 
where hi = ((0.2 - l/ 2 :)/(l + 1/z))^/^, /i 2 = ((0.2 -|- 
l/z)/(l-l/z))V3. 


> 



Figure 24. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,2 (black points) for the collection of func¬ 
tions [l,/ii,/i 2 ], where hi = ((0.2 — l/z)/(l -|- l/z))^/^, /12 = 
((0.2 + l/z)/(l-l/z))i/3. 
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Figure 25. Numerical distribution of zeros of type I Hermite- 
Pade polynomials < 5200,0 (blue points), < 5200,1 (red points), < 5200,2 
(black points), for the collection of functions [1, /i, / 2 ], where /i = 
log((0.4 - l/z)/(l + l/z)), /2 = log((0.4 + l/ 2 ;)/(l - l/z)). 



Figure 26. Numerical distribution of zeros of type I Hermite- 
Pade polynomials < 5200,0 (blue points) for the collection of func¬ 
tions [l,/i,/ 2 ], where /i = log((0.4 - l/z)/{I + l/z)), /2 = 
log((0.4 + l/z)/(l-l/z)). 
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Figure 27. Numerical distribution of zeros of type I Hermite- 
Pade polynomials Q 200 ,i (red points) for the collection of functions 
[1, /i, / 2 ], where /i = iog((0.4 - l/z)/{l + l/z)), /2 = log((0.4 + 
l/z)/{l-l/z)). 



Figure 28. Numerical distribution of zeros of type I Hermite- 
Pade polynomials ( 5200,2 (black points) for the collection of func¬ 
tions [l,/i,/ 2 ], where fi = log((0.4 - l/z)/(l -b 1/z)), /2 = 
log((0.4 + l/z)/(l-l/z)). 
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Figure 29. Numerical distribution of zeros of type I Hermite- 
Pade polynomials ^ 200,0 (blue points), ^ 200,1 (red points), ^ 200,2 
(black points), for the collection of functions [ 1 , 31 , 52 ], where 31 = 
((0.4 - l/z)/{l + l/z)Yl\ 32 = ((0.4 + l/z)/(l - l/z))i/2. 
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Figure 30. Numerical distribution of zeros of type I Hermite- 
Pade polynomials ^ 200,0 (blue points) for the collection of functions 
[ 1 , 31 , 32 ], where 31 = ((0.4 - l/ 0 )/(l + l/ 2 ))^/^ 32 = ((0.4 -f 
l/z)l{l-l/z)f/\ 
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Figure 31. Numerical distribution of zeros of type I Hermite- 
Pade polynomials P 200.1 (red points) for the collection of functions 
[ 1 , 51 , 52 ], where gi = ((0.4 - l/z)/(l + l/z))^/^ 52 = ((0.4 + 



Figure 32. Numerical distribution of zeros of type I Hermite- 
Pade polynomials P200,2 (black points) for the collection of func¬ 
tions [1,51,52], where 51 = ((0.4 - l/z)/{l + llz))^/"^, 52 = 

((0.4 + l/z)/(l-l/z))V2. 
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Figure 33. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,0 (blue points), 1 / 200,1 (red points), 1 / 200,2 
(black points), for the collection of functions [1, /ii, I 12 ]) where hi = 
((0.4 - l/z)/{l + /12 = ((0.4 + l/z)/(l - l/z))i/3. 



Figure 34. Numerical distribution of zeros of type I Hermite- 
Pade polynomials Q 200,0 (blue points) for the collection of func¬ 
tions [l,lii,/i 2 ], where hi = ((0.4— l/z)/(l -|- 1/z))^/^, I 12 = 
((0.4-bl/z)/(l-l/z))i/3. 
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Figure 35. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,1 (red points) for the collection of functions 
where hi = ((0.4 - l/ 2 ;)/(l + 1/z))^/^, /i 2 = ((0.4 + 
l/z)/(l-l/z))V3. 



Figure 36. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,2 (black points) for the collection of func¬ 
tions [l,/ii,/i 2 ], where hi = ((0.4— l/z)/(l -|- l/z))^/^, /12 = 
((0.4 + l/z)/(l-l/z))i/3. 
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Figure 37. Numerical distribution of zeros of type I Hermite- 
Pade polynomials <5200,0 (blue points), < 5200,1 (red points), < 5200,2 
(black points), for the collection of functions [1, /i, / 2 ], where /i = 
log((0.625 - l/z)/(l + 1/z)), /2 = log((0.625 + l/z)/(l - l/z)). 
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Figure 38. Numerical distribution of zeros of type I Hermite- 
Pade polynomials < 5200,0 (blue points) for the collection of func¬ 
tions [l,/i,/ 2 ], where ’/i = log((0.625 - l/z)/{I + l/z)), /a = 
log((0.625-bl/z)/(l-l/z)). 
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Figure 39. Numerical distribution of zeros of type I Hermite- 
Pade polynomials Q 200 ,i (red points) for the collection of functions 
[l,/i,/ 2 ], where/i = log((0.625-l/z)/(l-bl/z)),/z = log((0.625-t 
l/z)/(l-l/z)). 
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Figure 40. Numerical distribution of zeros of type I Hermite- 
Pade polynomials Q 200,2 (black points) for the collection of func¬ 
tions [l,/i,/ 2 ], where /i = log((0.625 - l/z)/{l + l/z)), /2 = 
log((0.625-bl/z)/(l-l/z)). 
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Figure 41. Numerical distribution of zeros of type I Hermite- 
Pade polynomials ^ 200,0 (blue points), ^ 200,1 (red points), ^ 200,2 
(black points), for the collection of functions [ 1 , 31 , 52 ], where 31 = 
((0.625 - l/z)/(l + l/z))i/2, 32 = ((0.625 + l/z)/(l - l/z))!/^. 
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Figure 42. Numerical distribution of zeros of type I Hermite- 
Pade polynomials ^ 200,0 (blue points) for the collection of functions 
[ 1 , 51 , 52 ], where 31 = ((0.625 - l/ 2 :)/(l + l/^))^/^ 32 = ((0.625 -t 
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Figure 43. Numerical distribution of zeros of type I Hermite- 
Pade polynomials P 200.1 (red points) for the collection of functions 
[1, 91 , 92 ]-, where gi = ((0.625 - l/ 2 :)/(l + 52 = ((0.625 + 

l/z)/{l-l/z)Y/\ 
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Figure 44. Numerical distribution of zeros of type I Hermite- 
Pade polynomials P 200,2 (black points) for the collection of func¬ 
tions [ 1 , 31 , 32 ], where 31 = ((0.625 - l/z)/(l -b l/z)Y/'^, 32 = 
((0.625+ l/ 0 )/(l-l/z))i/ 2 . 
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Figure 45. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,0 (blue points), 1 / 200,1 (red points), 1 / 200,2 
(black points), for the collection of functions [1, /ii, I 12 ]) where hi = 
((0.625 - l/z)/(l + l/z))i/3, h 2 = ((0.625 -t l/z)/{l - l/z)fl^. 
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Figure 46. Numerical distribution of zeros of type I Hermite- 
Pade polynomials 1 / 200,0 (blue points) for the collection of functions 
[1, hi, / 12 ], where hi = ((0.625 — l/z)/(l -I- 1/z))^/^, /12 = ((0.625 -I- 
l/z)/{l-l/z)Yl\ 
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Figure 47. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,1 (red points) for the collection of functions 
[1, /ii, / 12 ], where hi = ((0.625 — l/z)/(l -I- 1/z))^/^, /12 = ((0.625 -I- 
l/z)/(l-l/z))V3. 



Figure 48. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,2 (black points) for the collection of func¬ 
tions [l,/ii,/i 2 ], where hi = ((0.625 — l/ 2 :)/(l -I- 1/z))^/^, /12 = 
((0.625-bl/z)/(l-l/z))i/3. 
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Figure 49 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials <5200,0 (blue points), <5200,1 (red points), <5200,2 
(black points), for the collection of functions [ 1 , /i, /2], where /i = 
log(( 0.73 - + 1 /z)), /2 = log(( 0.73 + l/z)/(l - l/z)). 



Figure 50 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials <5200,0 (blue points) for the collection of func¬ 
tions [l,/i,/2], where /i = log(( 0.73 - l/z)/(l + l/z)), /a = 
log(( 0.73 + l/z)/(l-l/z)). 
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Figure 51. Numerical distribution of zeros of type I Hermite- 
Pade polynomials Q 200 ,i (red points) for the collection of functions 
[l,/i,/ 2 ], where /i = log((0.73 - l/z)/(l + 1/z)), /2 = log((0.73 + 



Figure 52. Numerical distribution of zeros of type I Hermite- 
Pade polynomials Q 200,2 (black points) for the collection of func¬ 
tions [l,/i,/ 2 ], where fi = log((0.73 - l/z)/(l + 1/z)), /2 = 
log((0.73 + l/z)/(l-l/z)). 
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Figure 53 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials ^200,0 (blue points), ^200,1 (red points), ^200,2 
(black points), for the collection of functions [1,31,52], where 31 = 
(( 0.73 - l/z)/(l + l/z))i/ 2 , 32 = (( 0.73 + l/z)/{l - l/z)Yl\ 
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Figure 54 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials P2oo,o (blue points) for the collection of functions 
[1,51,52], where 31 = (( 0.73 - l/z)/{l + l/z))^/^, 32 = (( 0.73 -t 

l/z)/(l-l/z))l/2. 
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Figure 55 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials P200.1 (red points) for the collection of functions 
[1,51,52], where gi = (( 0.73 - l/z)/(l + 32 = (( 0.73 + 

l/z)/(l-l/z))V 2 . 



Figure 56 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials P200,2 (black points) for the collection of func¬ 
tions [1,51,52], where gi = (( 0.73 - l/z)/{l -f g2 = 

(( 0 . 73 +l/z)/(l-l/z))i/ 2 . 
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Figure 57. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 300,0 (blue points), C/ 300,1 (red points), C/ 300,2 
(black points), for the collection of functions [1, /ii, / 12 ]) where hi = 
((0.73 - l/z)/{l + l/z))i/3, = ((0.73 + l/z)/{l - l/z))i/3. 



Figure 58. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/300,0 (blue points) for the collection of functions 
[1, /ii, / 12 ], where hi = ((0.73 - \/z)/{l + 1/z))^/^, /i 2 = ((0.73 + 
l/z)/{l-l/z)Y/\ 








ZERO DISTRIBUTION OF HERMITE-PADE POLYNOMIALS 


57 



Figure 59. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,1 (red points) for the collection of functions 
[1, /ii, / 12 ], where hi = ((0.73 - l/^)/(l + 1/z))^/^, /12 = ((0.73 + 
l/z)/(l-l/z))V3. 



Figure 60. Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,2 (black points) for the collection of func¬ 
tions [l,/ii,/i 2 ], where hi = ((0.73 — \/z)/{l + l/z))^/^, /12 = 

((0.73-tl/z)/(l-l/z))i/3. 
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Figure 61 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials < 5200,0 (blue points), <5200,1 (red points), <5200,2 
(black points), for the collection of functions [ 1 , /i, /2], where /i = 
log(( 0.8 - l/z)/{l + l/z)), /2 = log(( 0.8 -b l/z)/{l - l/z)). 
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Figure 62 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials <5200,0 (blue points) for the collection of func¬ 
tions [l,/i,/2], where /i = log(( 0.8 - l/z)/{I + l/z)), /2 = 
log(( 0.8 + l/z)/(l-l/z)). 
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Figure 63 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials Q200,i (red points) for the collection of functions 
[l,/i,/2], where /i = iog(( 0.8 - l/z)/{l + l/z)), /2 = log(( 0.8 -t 
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Figure 64 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials (5200,2 (black points) for the collection of func¬ 
tions [l,/i,/2], where fi = log(( 0.8 - l/z)/{l + l/z)), /2 = 
log(( 0.8 + l/z)/(l-l/z)). 
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Figure 65 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials ^200,0 (blue points), ^200,1 (red points), ^200,2 
(black points), for the collection of functions [1,31,52], where 31 = 
(( 0.8 - l/z)/{l + l/z)Yl\ 32 = (( 0.8 + l/z)/(l - l/z))i/ 2 . 
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Figure 66 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials ^200,0 (blue points) for the collection of functions 
[1,51,52], where 31 = (( 0.8 - llz)/{l + 1/2))^^^, 32 = (( 0.8 -t 
l/z)l{l-l/z)f/\ 
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Figure 67 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials P200.1 (red points) for the collection of functions 
[1,51,52], where gi = (( 0.8 - l/z)/(l + 52 = (( 0.8 + 
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Figure 68 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials P200,2 (black points) for the collection of func¬ 
tions [1,51,52], where 51 = ((0.8 - l/z)/{l + llz))^/"^, 52 = 

((0.8 + l/z)/(l-l/z))V2. 
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Figure 69 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/ 200,0 (blue points), 1 / 200,1 (red points), 1 / 200,2 
(black points), for the collection of functions [ 1 , /ii, I12]) where hi = 
(( 0.85 - l/z)/{l + l/z)fl^, h2 = (( 0.85 -b l/z)/{l - l/z))i/ 3 . 
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Figure 70 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials 1/200,0 (blue points) for the collection of functions 
[ 1 , 111 , /12], where hi = (( 0.85 - \/z)/{l + 1 /z))^/^, /i2 = (( 0.85 -b 
l/z)/{l-l/z)Y/\ 
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Figure 71 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials [/2oo,i (red points) for the collection of functions 
[ 1 , /ii, /12], where hi = (( 0.85 - l/^)/(l + 1 /z))^/^, /12 = (( 0.85 + 
l/z)/(l-l/z))V 3 . 



Figure 72 . Numerical distribution of zeros of type I Hermite- 
Pade polynomials C/200,2 (black points) for the collection of func¬ 
tions [l,/ii,/i2], where hi = (( 0.85 — \/z)/{l + l/z))^/^, /12 = 

(( 0 . 85 -tl/z)/(l-l/z))i/ 3 . 
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